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Abstract 

We calculated numerically the fractal dimension of the boundaries of the Fortuin- 
Kasteleyn clusters of the q-state Potts model for integer and non-integer values of q on 
the square lattice. In addition we calculated with high accuracy the fractal dimension of 
the boundary points of the same clusters on the square domain. Our calculation confirms 
that this curves can be described by SLE K . 
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1 Introduction 

Studying the critical interfaces in two dimensions is in the interest of both physicists and 
mathematicians because of at least two reasons: firstly because of application in the physical 
systems such as domain walls in statistical models [U [2] , iso- height lines in rough surfaces [3] , 
zero-vorticity lines of Navier-Stokes turbulence [4J and nodal lines of random wave functions 
[5]. Secondly these interfaces when they posses conformal symmetry can be studied by the 
exact methods of conformal field theory (CFT) [6| and rigorous methods of Schramm-Loewner 
evolution (SLE) [7]. 

Although physicists introduced many interesting results by considering the conformal in- 
variance of the interfaces there is no rigorous proof for most of them. It was just recently that 
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mathematicians were able to prove in some cases the conformal invariance of the statistical 
models at the critical point by using SLE techniques. Using SLE techniques Smirnov proved 
[8] the conformal invariance of percolation clusters which was conjectured long time ago by 
Cardy [9j and checked numerically by Langlands et al. |loJ^ . The proof was based on site 



percolation on the triangular lattice. There were lots of attempts to find a similar proof for 
other lattices but up to now the solution is out of reach [11]. In 2007 again Smirnov proved 
the conformal invariance of Ising clusters in the Fortuin-Kasteleyn (FK) representation on 
the square lattice [12J. Finally recently Chelkak and Smirnov proved the conformal invariance 
of spin clusters of Ising model on the generic iso-radial graph [13] . Apart from the above 
list of rigorous proofs of well-known statistical models there are just three more proofs for 
conformal invariance of stochastic curves, loop erased random walk [14] , harmonic explorer 
[15] and the contour lines of Gaussian free field theory [16] . 

Potts model is one of the building blocks of statistical mechanics and it was studied from 
many different points of view. In two-dimensions at the self dual points it is an exactly 
solvable model and the critical properties were studied in detail, see [17] . It was also studied 
from CFT point of view [18J. The critical domain walls in Potts model were also studied 
in detail by using Coulomb gas techniques and lots of conjectures came out from heuristic 
arguments [El [20]. Some numerical calculation were done to support the conjectures for 
the Potts models with integer q, see (2TJ [22j [231 12H ES] . They mostly calculated the fractal 
dimension of domain walls, cluster masses, red bonds, fjords and singly connected bonds, by 
simulating them on the square lattice in the bulk of the rectangular domain. Although none 
of them separately guarantees that those curves are SLE and so conformal since there are 
some relations between the above quantities and k of SLE for conformal curves [23] . one can 
get convinced that they are related to SLE after checking at least the validity of two of them. 
In other words if two of the above quantities satisfy specific relation between each other then 
the curves are conformal. The above argument still is not enough to rule out for example the 
possibility of having SLE(rc, p), see [26]. To check that a process is SLE(k) or SLE(k, p) one 
needs to check directly the properties of the drift of SLE or to check some boundary fractal 
properties. It is difficult to get some results with high accuracy by working directly with SLE 
equation [27]. However, it is easier to calculate the boundary fractal dimension of the curves 
in simulation. 



1 These authors checked numerically the invariance under Mcebius transformation and not invariance under 
the whole conformal maps. 
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By the above motivations in this article we will give some strong non-direct simulation 
supports for the conformal invariance of the boundary of FK clusters in the Potts model for 
integer and non-integer values of q. To the best of our knowledge for the non-integer case 
there is no simulation or rigorous argument for having conformal domain walls. To rule out 
the possibility of having SLE(k, p) we will calculate the boundary fractal dimension of the 
curves. 



2 Definition and Theory 

The definition of g-state Potts model on the square lattice in the arbitrary domain is as 
follows: associate a spin variable Si 6 {0, 1, q — 1} at each site then the partition function 
is 

Z = £ ^ = E II (! + s i))> 

Si s-i <i,j> 

where u = e J — 1. By expanding the product one can obtain the FK representation as follows 

£ = 5>V, (2.2) 

G 

where G is any subgraph of the original domain, consisting of all the sites and some bonds 
placed arbitrarily on the lattice edges, b is the number of bonds in G, and c is the number of 
clusters of connected sites into which the bonds partition the lattice. Although the original 
formulation of the model requires q to be a positive integer, the FK representation allows one 
to interpret q as taking arbitrary real values. The model is known to have a critical point at 
the self-dual value of u c = ^fq for < q < 4. The partition function can be also expressed 
as that of a gas of fully packed loops on the medial lattice. At the critical point FK clusters 
are equivalent to counting each loop on the medial lattice with a fugacity u c , see |28j . It was 
conjectured in [29] that these loops can be described by SLE with the following equation 

4tt 

« = (2-3) 

cos-i(-^) 

where 4 < k < 8. Since the fractal dimension of SLE is 

<*/ = 1 + f> ( 2 - 4 ) 
(see |30|). then one can find the following formula for the fractal dimension of loops in the 
q-state Potts model 
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This is a well known formula in physics literature from the Coulomb gas arguments [20] but 
it has not been checked numerically for non-integer values of q. In the next section we will 
check the validity of the above equation for some values of q in the rectangular domain. The 
popular version of SLE is defined on the upper half plane, a stochastic curve emerges from 
origin and goes to infinity by keeping the right-left symmetry. Since SLE is a conformal 
curve, mapping the upper half plane to rectangle does not change the fractal properties of 
the curve; for the non-conformal curves this argument is not true. With the above argument 
one should be careful that checking the equation (|2.5h on the rectangular domain does not 
say anything about the conformality of the curve. As we argued in the introduction checking 
at least one more quantity is necessary for stronger numerical argument. One interesting 
quantity is the fractal dimension of the boundary points of the fractal curve. It was proved 
in |31| that the fractal dimension of the SLE points on the real line is 



then one can find the following formula for the fractal dimension of boundary points for the 
q-state Potts model loop 



Checking the above equation in the rectangular domain is a stronger numerical check for 
conformality of the stochastic curves. The above equation is also important to support that 
these curves are related to SLE(k) and not to SLE(k, p). The latter one is also a conformal 
curve and has similar bulk properties but the fractal dimension of the boundary points is 
different [32]. The most famous statistical model related to SLE(k,/?) is the contour lines of 
Gaussian free field theory [16] . 

We will close this section by a short comment on the relation of the bulk and boundary 
fractal dimensions to the conformal field theory operators. From CFT and Coulomb gas 
arguments it is well known that ^>o,f 1S related to the point with n attached curves. Since the 
bulk fractal dimension is related to n = 2 then the responsible bulk field is </>o,i with weight 
2/io,i = 1 — f • Then the fractal dimension is d = 2 — 2/10,1. To produce a boundary point 
with n attached curves one needs to plug the operator <^i, n +i at the corresponding point. 
Then the fractal dimension of boundary points is = 1 — /ii 5 3, where h\ 3 is the weight of 
the operator <f>i$. 




(2.6) 



K 





TT 



4 



Combining equations (|2.4p and (|2.6p one obtains the relation 

(d / -l)(2-d 6 ) = l , (2-8) 

where any specific reference to the q-state Potts model has disappeared, therefore one could 
suspect that such a relation between the Hausdorff dimensions df and df, could be true for a 
wider class of fractal curves. 

In the next section we will give some extensive numerical simulations to check the validity 
of the equations (|23]) . (I2T7D and ([22]). 



3 Simulations 
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Figure 1: Setting of our numerical simulations: in a square box of size L x L we measured 
the sum of the perimeters (dashed line)of all the FK clusters (thick lines) touching a fixed 
segment of the boundary of length I = L/A placed symmetrically at the centre of one side. 
We also evaluated the fractal dimension of the boundary by counting the number of the 
boundary sites (black dots) touched by these clusters. 



We simulated three different systems in a square lattice: a bond percolation system at the 
threshold, corresponding to a q = 1 Potts model with k = 6, a critical q = 2 Potts associated 
to k = j and finally a critical q = 2cos 2 (^) = "^g 3 Potts model which according to (|2.3p 
corresponds to k = 5. In all cases the system was enclosed in a square box of size L x L. 



We measured the sum of the perimeters of all the FK clusters touching a segment of length 
£ = L/4 placed symmetrically on the centre of a side of the box (see Figured]). 

In the case of percolation the boundary conditions were chosen free along the whole 
perimeter B of the square box. We adopted an epidemic type of algorithm: to begin the 
process, all the sites of the lattice are set to the "unvisited" state and bonds to the "unde- 
termined" state. Then we pick one of the unvisited sites i of the segment £ and we check 
the three bonds connected to i; the undetermined bonds are set "occupied" with probability 
P = Pc = \ and "empty" otherwise. For bonds that are occupied, we check the adjacent site; 
if that site is unvisited we label it as visited and put it in a list for further checking. After 
finishing checking all bonds connected to a site, one considers the next site on the list, contin- 
uing this process until the list is empty. One repeats this process for each unvisited site of the 
segment £. In this way we generate only the clusters emanated by £. This algorithm allows 
to consider boxes of large size with little computational effort. In each sample generated this 
way we measured the sum of the perimeters of these clusters as well as the number of sites 
of the boundary reached by these clusters. We generated two sets of data; in the first set 
we considered lattices of size L x L with L = 80, 120, . . . , 720 with 1.2 10 7 samples for each 
L. The second set was composed of larger boxes (L = 1360, 1600, 1800, 2000, 2400, . . . ,4000) 
with about 6 10 6 samples for each L. Figure [2] shows the results for a plot of the total 
perimeter of the FK clusters versus L in a logarithmic scale 

In the cases of Potts models with q > 1 we performed two different kinds of Monte 
Carlo simulations. For the estimate of the fractal dimension df of total perimeter of the FK 
clusters emanated by the segment £ we chose fixed boundary conditions along this segment, 
say Si = Vi G £, and Sj^OViG B \ £, while in estimating the Hausdorff dimension 
db of the boundary we chose always fixed boundary conditions along £, but free boundary 
conditions along the rest B \ £ of the square box. While these boundary conditions have 
an obvious implementation in the q = 2 case where we adopted a standard Swendsen-Wang 
algorithm [33], they need further specification in the case of non integer q, where the spin 
variables Sj are ill-defined. Among the proposed algorithms for non integer q [331 ESI ES] , the 
most efficient one for the q > 1 case is that of Chayes and Machta, which we implemented in 
the following form, starting from an arbitrary configuration of clusters of connected sites 

1. set all clusters to the "undetermined" state; 

2. assign to each undetermined cluster a white color with probability | and a non- white 
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Figure 2: Plot of the total perimeters of the clusters emanated by the segment I as a function 
of the size L of the percolation process. The data lie on a straight line in this logarithmic 
plot. The slope yields an estimate of the fractal dimension df = 1.74990 ±0.00002 in perfect 
agreement with the expected exact value df = 1 + |. 



1000 




10000 



L 



Figure 3: Plot of the mean number of boundary sites attained by the FK clusters in the 
three Potts models studied. The three continuous lines are the fits of the data to aL db . The 
size of the systems with q > 1 is much smaller than in the case of percolation causing larger 
statistical errors in the determination of df,, as Table 1 shows. 
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color otherwise; 

3. set each bond between white sites "occupied" with probability an d "empty" oth- 
erwise. In this way a new configuration of clusters of connected sites is generated; 

4. return to step 1. 

Table 1: Estimated fractal dimension df and df, in the critical Potts models with q = 1, q = 2 
and q = v/ ^ +3 compared with the expected value in SLE. 
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df 


expected df 


d b 


expected df, 
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1.74990 ± 0.00002 


7 
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0.6656 ± 0.0002 


2 
3 
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1.659 ±0.004 


5 
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0.505 ± 0.003 


1 

2 


v^+3 
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1.625 ±0.003 


13 
8 


0.41 ±0.01 
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Note that this algorithm uses only two kinds of sites - white and non-white - but they are 
not treated equivalently. In our numerical estimates of df we chose white sites in the segment 
£ and non- white sites in the rest of the boundary B\£, while in the case of the determination 
of the fractal dimension d b the boundary no restriction was made on the sites of B \ I. 

In order to extract our numerical estimates of the fractal dimension we fitted our data to 
the power law a L d using a and d as fitting parameters. Since this is an asymptotic expression, 
valid when the size L of the box is much larger than the lattice spacing, we fitted the data 
to the power law by progressively discarding the short distance points until the x 2 test gave 
a good value. The estimated values are reported in Table 1 and in Figured! 

4 Conclusions 

In this paper we calculated numerically the fractal dimension of the boundaries of the FK 
clusters of the Potts model, for the integer and non-integer values of q, on the square lattice. 
In addition we found precisely the fractal dimension of the boundary touching points of the 
FK clusters on the square lattice. Since the bulk fractal dimension and the boundary fractal 
dimension of the domain walls of FK clusters of the Potts model have the same relations 
that we expect from SLE we believe that our method gives convincing numerical support for 
the conformal invariance of the boundaries of the FK clusters. The direct method to see the 
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Figure 4: Plot of the estimated fractal dimensions in the plane {db, dt). The continuous curve 
represents Eq. (|2.8|) . 

relation between the critical curves and SLE is calculating the drift of the Loewner equation 
and checking its equality, in the distribution sense, with the Brownian motion. The numerical 
methods in this direction is not that much efficient and one can not get very precise numbers. 
For this reason checking at least two fractal properties of the curves can give more efficient 
numerical support for the conformal invariance of the model. 

Acknowledgment: We are indebted to M. Caselle and S. Lottini for fruitful discussions. 

Not added: After completion of the present work, the paper |25j has appeared. It also calcu- 
lates the fractal dimension of the boundary of FK clusters with the comparable accuracy for 
integer q's. 

References 

[1] M. Bauer and D. Bernard, Phys. Rep. 432 (2006) 115 [ }ArXiv:math-ph/0602049] , J. 
Cardy, Annals Phys., 318 (2005) 81 |ArXiv:cond-mat /0503313] . 

[2] C. Amoruso, A. K. Hartmann, M. B. Hastings, M. A. Moore, Phys. Rev. Lett. 97 (2006) 
267202 [ArXiv:cond-mat/06017Tl |; J. L. Jacobsen, P. Le Doussal, M. Picco, R. San- 




9 



tachiara, K. J. Wiese, Phys. Rev. Lett. 102 (2009) 070601 [arXiv:0809.3985] : M. Picco, 
R. Santachiara, A. Sicilia, J. Stat. Mech.(2009) P0401 |arXiv:0812.3526] : M. Caselle, S. 
Lottini, M. A. Raiabpour. larXiv:0907.5094[ 

[3] A. A. Saberi, M. A. Rajabpour, S. Rouhani, Phys. Rev. Lett. 100 (2008) 044504 
[arXiv:0712.2984j . 

[4] D. Bernard, G. Boffetta, A. Celani, G. Falkovich, Nature Phys. 2 (2006) 124 
[ArXiv:nlin/0 602017|. 

[5] J. P. Keating, J. Marklof, I. G. Williams, Phys. Rev. Lett, 97, 034101 (2006) 
[ArXiv:nlin/06030 68| . 

[6] A. A. Belavin, A. M. Polyakov, A. B. Zamolodchikov, Nucl. Phys. B 241 (1984) 333. 

[7] O. Schramm, Israel J. Math. 118 (2000) 221 [ArXiv:math/9904022| . 

[8] S. Smirnov, C.R. Acad. Sci. Paris, 333 (2001) 239 garXiv: 0909 .4499] . 

[9] J. L. Cardy, J. Phys. A 25 (1992) L201 [ArXiv:hep-th/9111026| . 

[10] R. Lang.lands, P. Pouliot, and Y. Saint-Aubin, Bull. Amer. Math. Soc. (N.S.) 30 (1994) 
1 |ArXiv:math/9401222| . 

[11] V. Beffara. la7Xiv:0708.3908l . 

[12] S. Smirnov, arXiv:0708.00"39l . 

[13] D. Chelkak, S. Smirnov, arXiv:0910.2045l . 

[14] G. F. Lawler, O. Schramm, W. Werner, Annals Probab., 32 (2004)939 
[ArXiv:math/01 12234] . 

[15] O. Schramm, S. Sheffield, Annals Probab., 33 (2005) 2127 [ArXiv:math/0310210| . 

[16] O. Schramm, S. Sheffield, Acta Mathematica, 202 (2009) 21 [ArXiv:math /0605337] . 

[17] R. J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic Press, New 
York, 1982). 

[18] P. Di Francesco, P. Mathieu and D. Senechal, Conformal Field Theory (Springer, 1997). 
[19] H. Saleur, B. Duplantier, Phys. Rev. Lett, 58 (1987) 2325. 

10 



[20] A. Coniglio, Phys. Rev. Lett. 62 (1989) 3054. 

[21] J. Asikainen, A. Aharony, B. B. Mandelbrot, E. M. Rauch, J.-P. Hovi, Eur. Phys. J. B 



34 (2003) 479 ArXiv:cond-mat/0212216 ]. 



[22] B. Wieland, D. B. Wilson, Phys. Rev. E 68 (2003) 056101 |arXiv:1002.3 220j. 



[23] W. Janke, A. M. J. Schakel, Nucl. Phys. B, 700 (2004) 385 |ArXiv:cond-mat/0311624|. 

[24] Y. Saint-Aubin, P. A. Pearce, J. Rasmussen, J. Stat. Mech. (2009) P02028 
[arXiv:0809.4806j . 

[25] D. A. Adams, L. M. Sander, R. M. Ziff, larXiv:1001.0055l 

[26] G. Lawler, O. Schramm, W. Werner, J. Am. Math. Soc. 16 (2003) 917 



[ArXiv:m ath /0209343] . 



[27] T. Kennedy, J. Stat. Phys. 128 (2007) 1125 [ArXiv:mat h /0508002] ] . 
[28] B. Nienhuis, J. Stat. Phys, 34 (1984) 781. 



[29] S. Rohde, O. Schramm, Annals Math. 161 (2005) 879 [ArXiv:m ath/0106036|. 

[30] V. Beffara, Annals of Probability, 36 (2008) 142 [ArXiv:math/0211322| . 

[31] O. Schramm, W. Zhou. larXiv:0711.3350l : T. Alberts, S. Sheffield. llarXiv:071 1.40701 . 

[32] J. Dubail, J. L. Jacobsen, H. Saleur, Phys. Rev. Lett, 103 (2009) 145701 
[arXiv:0909.2949| 

[33] R.H. Swendsen and J.-S. Wang, Phys. Rev. Lett. 58 (1987) 86. 

[34] M. Sweeny, Phys. Rev. B 27 (1993) 4445. 

[35] L.Chayes and J. Machta, Physica A 254 (1998) 477. 

[36] F. Gliozzi, Phys. Rev. E 66 (2002) 016115 [ArXiv:cond-mat /0201285] . 



11 



